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Abstract 

Parton recombination is reconsidered in perturbation theory without using the 
AGK cutting rules in the leading order of the recombination. We use time-ordered 
perturbation theory to sum the cut diagrams, which are neglected in the GLR 
evolution equation. We present a set of new evolution equations including parton 
recombination. 
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1 Introduction 



Parton recombination as a new higher twist phenomenon was first discussed in the 
QCD evolution process by Gribov, Levin, Ryskin [1] and Mueller, Qiu [2] in their pio- 
neering works. This evolution equation is called the GLR equation. 

An interesting effect of parton recombination is screening or shadowing. In the case 
of higher number densities of partons, for example in the small x region, the gluons can 
overlap spatially and annihilate. Therefore, one expects that the growth of the gluon 
density with Q 2 will be suppressed by gluon recombination. These suppression factors 
from the negative contributions due to gluon recombination are calculated in the GLR 
equation using the AGK (Abramovsky, Gribov, Kancheli) cutting rules [3]. Assuming 
that the AGK cutting rules are valid in deep inelastic scattering (DIS) in the small x 
region, one finds that the relative weights of cuts through two, one and zero ladders are 
2 : —4 : 1, as illustrated in fig. 1. Due to these quantitative predictions of the suppression 
of parton number densities at small x, the GLR equation was extensively used to explore 
the structure of the nucleon and new perturbative QCD (PQCD) effects in the past years. 

However, the applications of AGK cutting rules in the GLR equation have some draw- 
backs. For example, the cut lines break the correlation between the recombining partons 
according to the AGK cutting rules in the GLR equation (fig. 1). As we will show in 
this work, the correlation among the initial partons in the QCD recombination equation 
should be preserved. On the other hand, two-to-two parton processes may be associ- 
ated with IR-divergences just as are the one-to-two processes in the Dokshitzer-Gribove- 
Lipatov-Altarelli-Parisi (DGLAP) evolution equation [4-6]. Although the AGK cutting 
rules provide simple relations between the cross sections of hadron-hadron interactions 
for different types of reggeon cuts, however, the sum of cut graphs according to the AGK 
cutting rules cannot cancel IR-divergences. The reason is that the difference of the con- 
tributions between the positive graph and the negative graph is only a weight, according 
to the AGK cutting rules. As we know that the virtual diagrams are necessary for cancel- 
lation of IR-divergences in the DGLAP equation. We will show in this work that above 
mentioned IR-divergences in two-to-two processes also can be canceled by the sum of 
virtual diagrams, which are neglected in the GLR equation. 

In this work, we reconsider parton recombination in the QCD evolution equation with- 
out the AGK cutting rules. To this end, we first point out that a new scale (the recom- 
bination scale) exists in the parton recombination processes. We will give a definition of 
the recombination order of the process. Then, we propose the bare probe-vertex approx- 
imation. We find that several more diagrams, which are neglected in the GLR equation, 
should be included in the QCD equation of the parton fusion. We try to find a simple way 
to calculate those cut diagrams in time-ordered perturbation theory (TOPT). Through a 
new derivation of the DGLAP evolution equation, we present simple connections among 
the different cut diagrams. We shall show that both the shadowing-antishadowing and 
momentum conservation are the natural results of the theory. As an interesting result, 
our new equation has different structure from the GLR equation. 

The outline of the paper is the following. In section 2 we give some definitions related 
to the parton model. In section 3 the sum of the cut diagrams at the bare probe-vertex 
approximation is proposed. In section 4 we give a new derivation of the DGLAP equation. 
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Through this example, we try to show the connections among the relative cut graphs at 
the leading recombination order. The new evolution equations incorporating the par- 
ton recombination are derived in sections 5-7. Section 8 contains the discussions and 
concluding remarks. 

2 Definitions 

According to PQCD, a parton can always independently split into two partons. How- 
ever, except in the PQCD dynamics, the recombination of partons depends on the overlap 
probability of their wave functions. Therefore, we need a new physical quality to charac- 
terize parton recombination. 

For example, consider an amplitude with two initial partons. According to dimensional 
analysis, the hadronic part of the amplitude should contain a factor ~ l/R, where R has 
dimension of length. Now we incorporate the factor 1/Q arising from the partonic part 
of the amplitude to form a dimensionless quality 1/(RQ). We call this the recombination 
factor. For example, the recombination factor of the process with two initial correlated 
partons is 1/(RQ) 2 . One can schematically think of 1/(RQ) 2 as the overlap probability 
of two partons, where 1/Q is the scale of a parton at momentum transfer Q 2 and R is 
the maximum correlation length of two partons. Usually, R is regarded as the scale of 
target or the scale of the "hot spots" , if they exist in the proton. This definition can be 
generalized to the case of the amplitudes including m-initial partons; the recombination 
factor in this case is l/(RQ) m ~ l if the fusing partons are paired. 

In this paper, we consider only the recombination processes at the leading order level, 
that is, at 1/(RQ) 2 and a 2 . Therefore, we choose to study the following basic amplitudes 
as shown in fig. 2: (a) M^^ k , vx , (b) M^7^ k , vx and (c) M^7^ k , llx , with the recombi- 
nation factors 1/(RQ)°, lf(RQ) 1 and 1/(RQ) 2 , respectively. In fig. 2 we have omitted 
the distinction of the parton flavors; the dark circles indicate QCD interactions among 
the correlating partons. 

We see that the amplitudes involving parton recombination contain the double scales: 
1/(RQ) and a s . We shall perform the calculations at a given order- l/(RQ) m and order- 
{a s ) n in two steps. First, we calculate the process at the order of 1/(RQ) 2 . In this case, 
the dark circles in fig. 2 are regarded as the elemental sub-processes. The reason is that 
the decomposition of the circle part will break the parton correlation and reduce the order 
of the recombination. The second step is to calculate the sub-processes at order-a;; in 
PQCD. 

We will use TOPT in this work. Usually, TOPT is equivalent to the standard covariant 
perturbation theory [7]. We shall call this TOPT as the normal TOPT (NTOPT), where 
the time lines divide every basic vertex along the time-order. In TOPT, the internal 
lines and the virtual particles are expressed by external lines and effective real particles 
respectively. Therefore, TOPT can also be used to describe amplitudes involving complex 
vertices, where the part between two neighboring time lines can contain a complex vertex. 
We define such a TOPT as an anomalous TOPT (ATOPT). Obviously, ATOPT is not 
equivalent to NTOPT: they have different energy deficits. On the other hand, there 
are energy-momentum correlation between two neighboring complex vertices in ATOPT; 
therefore, the vertex in ATOPT is not really factorized. 
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We take the physical axial gauge, where the light-like vector n fixes the gauge as 
n • A = 0, A being the gluon field. The parton number densities are defined within the 
parton model description of the photon nucleon DIS (fig. 3a) as 

da(-f*p -> k'X) = J dx 1 q{x 1 )da{^k -> k'), (2.1) 

where q(x\)dxi is the number of quarks carrying momentum fraction between x\ and 
x\ + dxi, where x\ = k ■ n/p ■ n. Formula (2.1) means that the interaction of a virtual 
photon with proton can be factorized as the soft part q(x\) and hard part da(j*k — > k'). 
According to the parton model, 

da(-f*k -> k') = C q 5( Xl - x B ), (2.2) 

where xb = Q 2 /2p.q and C q is the coefficient depending on x B and Q 2 . The quark 
density can be defined as 

q( XB ) = _L da^p^k'X). (2.3) 

Q 

On the other hand, using TOPT in the cut graph 3b, we have 

da(-f*p -> k'X) 

E k lA , ,2r 1 i2r 1 12" d% 



- -rr\M p ^ kX \ 2 [ ] 2 [ ] 2 J] 



-Ep E p — E k — E x 2E k x (2n) 3 2E x 

'M rfc ^| 2 (27r)V(p 7 + fc - k') —^— . (2.4) 



8^' v y (2vr)3 2 ^ 

Comparing eqs. (2.1) with (2.4), we get the definitions of the quark number density 

q{x\)dx\ 

= fj Mp - kXn E p - E k - E f [ ^ ]2 J ^e~ x ( 2 - 5 ) 
and the bare probe-parton vertex 

^rda{Yk -> fc') 
|M 7 ^| 2 (27r) 4 ,5 4 (j9 7 + A;-A;') 



C q 8E k E^ »k , v^ 7 ^(27r) 3 2^ 

= 5(zi-zb), (2.6) 

in the TOPT-form. 
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3 Bare probe-vertex approximation 



As we know, emission or absorption of quanta with zero-momentum may associate with 
the infrared (IR) divergence. However, the singular terms provide the leading contribu- 
tions to the DIS processes. Therefore, a correct theory is IR-safe, where IR-divergences 
are canceled, while the leading contributions are retained. One way can to attain above 
two goals is to sum over cut diagrams belonging to the same time-ordered uncut graph, 
since these graphs have similar singular structure but may come up with opposite signs. 

Deep inelastic scattering structure functions are the imaginary parts of the amplitudes 
for the forward 'Compton' scattering of the target with a probe. Using the time-ordered 
perturbative expansion of the statement of the unitarity of the S-matrix, one can prove 
that the structure functions are associated with the sum of cut diagrams. These different 
cut graphs represent various possible sub-partonic processes due to the unitarity of the 
perturbative S-matrix [7]. Therefore, the sum of cut graphs is necessary not only for 
infrared safety, but also for collecting the leading contributions and restoring the unitarity. 

The interesting and important question is, what are the minimum cut diagrams that 
must be summed for IR-safe calculation of an inclusive DIS process at a given order (for 
example, order a 2 s /(RQ) 2 in this work) ? 

To answer this question, let us consider a general inclusive DIS process on target N. 
One can choose the cut diagrams according to following program: G(N) stands for the 
time-ordered uncut diagram of the target N without probe vertices. We sum over possible 
cut diagrams of G(N): 

£G 7 (A0 = {£V*7}/, (3-1) 

7 7 

where G 7 is the cut diagram with cut line 7; L 7 and i? 7 are the sub-graphs on the 
left and right of the cut line; the subscript "I" means that we only consider following 
cut graphs which have the same observed quantities (that is, the same structure of the 
intermediate state) and which keep the original correlation among initial partons in G(N). 

We use the probe to observe the parton distributions inside the target in DIS. Of 
course, we cannot control the probing positions. In principle, in- and out-probe lines 
can be attached to the left- and right-hand of the cut line in all possible ways. Let 
G@ (probe + N) stand for the cut diagrams of the probe-target system, where (3 labels the 
probe-parton vertices. Thus, we shall sum over 

Y / J2 G ^(P robe + N )- ( 3 - 2 ) 

p 7 

Obviously, the sum (3.2) is much larger than J2- y G 1 (N). Now we try to find some 
approximation in the sum (3.2). As we know that the leading logarithmic approximation 
(LLA) is a good approximation for IR-safeness in the DGLAP equation at order-o; s . In 
this approximation, some of the renormalization effects are neglected in the physical gauge 
and the probe- vertex retains the bare- vertex form as in (2.6). In this case, 

G^probe + N) = G 1 {N)8{x 1 - x B ), (3.3) 
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Thus, the contributions from the nonlocal interactions of probe with partons are ne- 
glected at the leading approximation. We need only to sum part of cut graphs, in which 
the bare-probe vertex 5(x 7 — x B ) connects with the cut line, that is, 

7 
7 

= J2i L ~f R ~f}i 5 ( x ~f - x b), (3-4) 

7 

in the DIS processes with parton recombination. 

We call (3.4) as the bare probe- vertex approximation. We find that this approximation 
is also a satisfactory approximation in the DIS processes with parton recombination. In 
fact, our interest is that the modifications of the parton recombination to the DGLAP 
equation, which has the probability explanation at the LLA approximation. We shall 
show that the bare probe-vertex approximation is necessary for keeping the probability 
picture of the new evolution equation. 



4 Rederivation of the DGLAP equation 

We know that several methods can be used to derive the DGLAP evolution equations, 
however, the following new method illustrates more clearly the simple relations among 
the cut diagrams in the sum (3.4). For simplicity, we only consider the non-singlet case. 
According to (3.4) we compute fig. 4a with figs. 4b and 4c. Figures 4b and 4c seem to 
change the observed quantity dlnl"j_, where l± is the transverse momentum of the final 
state partons, (see fig. 2a) whose momenta are parametrized as 



I = [xip, 0, xxp] , 



k = 



I' 



x 2 p + 
x 3 p + 



V 



2x 2 p 
I 2 



,l±,x 2 p 



2x 3 p 



, -l±,x 3 p 



(4.1) 



However, we will see that the contribution of an intermediate state can be replaced 
by d 3 l', which is from the contribution of the loop and contributes d\nl 2 ± . Therefore, all 
the processes of fig. 4 have the same intermediate state structure. 

We proceed along the lines of ref . [8] . The change of the valence-quark-number density 
caused by gluon radiation can be written as (see fig. 4a) 



dq(x B ) 



a 



■da(<y*p -> k'X) 



E 



E n 



Ip-*IX\ 



d 3 k, 



E n — E\ 



— 1 2 [— i 2 n 

E x l [ 2E t l V (2tt)32£ x 
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ff( 7 */->7*0, (4-2) 

where the cross section is factorized to the soft part and the hard part H{^*1 — > 7*/) 
according to the factorization theory [9,10]. Using TOPT we obtain the hard part 

H ( 7 */ - 7*/) 
1 ' ■\M rk ^ k/ \ 2 (27t) 4 5\p, + k-k')- ,n ' 



C q 8E k E^ > ^7 J (2nf2E k , 

Assuming that /j_ ~ Q 2 , we have 

cig(xB,Q 2 ) 
dlnQ 2 

= J 9 ( Xu Q^P^ *,)*(*! - - 2 - X 8 )^ a - X^MXA 

= I q{x u Q 2 )Pn-)—, (4.4) 
7 xi Xi 



where 



qM * )dxi = I.M^ri^^-J^]^]^^-, (4.5) 



as same as eq.(2.5) and 



I I rl 3 k' 

In eq.(4.4) we inserted 



y —x 2 — x 3 )dx 2 = 1. (4.7) 



We define 



J/2 771 1 1 j3// 

^.^^f = ||^|'[^-A_]. [5 Lp^., (4.8) 
as the parton splitting function for the non-singlet part. 

Now let us consider figs. 4b and 4c. As for the the real diagram, the contributions of 
fig. 4b to the change of the valence-quark-number density are 

dq{x B ) = \~-do{l*V -> k'X) 

-^n P _:,_ g f[^n^ d3kx 



2E P < p lA 1 V E V - Ei - Ex 1 L 2^ J A / (2vr) 3 2£ x 
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Htfl^^l), (4.9) 

where the factor of \ was explained as the effect of the renormalization in ref. [11]: 
The virtual part of fig.4b corresponds to the renormalization of a parton propagator. Only 
half of the probe-vertex connects with this parton line. This is equivalent to multiplying 
the virtual process by an extra factor of |. 

The hard part in eq.(4.9) is 

Htfl - 7*0 
' \M lH ^ k ,\ 2 (2vr)V(p 7 + /-A;') " 



C q mEj 71 ^ v ' J (2irf2E k , 

d* 1 ' 1 1 1 1 1 

rMi^ky——- —M kl ,_+ r 



(2tt) 3 2E k 2E v E l -E k -E v E k + E v — E\ 2E, 

1 1 rl 3 h' 

\M^ kl \\2nY5\ Pl + l-k') ' 



C q mE^ 7 '^' v 1 ^ '(2ixf2E k , 

^\ M ^ kl ,\ 2 [ Ft—] 2 , f/ ■ (4.10) 

El \ [ El -E k -E v ] [ 2E k ] {2ixf2E v 1 ; 

Therefore we have 

dq(x B ,Q 2 ) 
dlnQ 2 

= - J ^q(x 1 ,Q 2 )P$ q (xi,x 2 ,x 3 )5(x 1 -x 2 - x 3 )8(xi - x B )dx 1 dx 2 dx 3 

= - J ^q(xi, Q 2 )P 2 q9 (z)dz5(x 1 - x B )dxi 

= -\q(x B ,Q 2 )Jpr(z)dz, (4.11) 
in which we have used eqs.(2.5),(2.6) and (4.8). 

The contributions of fig. 4c is same as one of fig. 4b. Thus, the total contributions of 
the real- and virtual-diagrams are 

= j q(xM ^ _ M) j P?(z)dz , (4 , 2) 

in which 

P?{z) = ^C 2 {R) l -±^. (4.13) 

The two terms of the right-hand side in eq.(4.12) have a simple interpretation: the 
positive contribution arises from the splitting of higher momentum quarks, while the 
negative contribution results in the loss of the number of quarks due to its gluon radiation. 

The result (4.12) is the same as the probability form of the DGLAP equation for the 
nonsinglet part in ref. [11]. However, the new derivation clearly shows the following inter- 
esting properties in the inclusive DIS processes: The contributions of the cut diagrams, 
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which belong to a same time-ordered uncut graph in the sum (3.4), have an identical 
integral kernel (it is the parton splitting function in eq.(4.12)). This is a reason we use 
the TOPT form to perform our calculations. We shall examine this connection further in 
the parton recombination processes. 



5 Leading recombination approximation 



So far we considered processes contributing to the usual DGLAP equation. We go on 
n to include parton recomt 
leading order come from tl 
diagrams according to (3.4). 



now to include parton recombination processes. The recombination processes contributing 
at leading order come from the terms, \M^I^ krvx \ 2 , 2M^7^ kq/x [M^ S l kll/x }* , an d the cut 



In this section we regard the partons as scalar particles (i.e., the 3 model). The 
results can easily be generalized to the case of QCD partons and will be done later in 
section 7. We consider the process of fig. 5, where figs. 5c-f are virtual. The contribution 
of the real diagram fig. 5a is 

d<P{x B ) = ^rda(-f*p -> k'X) 



^E Pl +E P2 ^E p[ +E p , 

— M p ^ PlP2 x [M p ^ p > iP ' 2 x\ 



, 2 1 1 I 1 j j d 3 k x 



Ep — E Pl — E P2 — Ex 2E Pl 2E P2 2E p > i 2E P > 2 x (2n) 3 2Ex 

H^Pm^YM), (5.1) 

where is defined by 

da(-f*k -> k') = C (j} 5(x 1 - x B ), (5.2) 
for the scalar parton; the hard part is given by, 

#(7*PiP2 -> 7V1P2) 
= (b 2 ^r^r\M rk ^\^m P , + k - k') (H '' 



R ; C^E k E^ k—>k j v^t '(2Tif2E kl 

:M pip2 _ >W /[Mp/ i p/_ >fc //]* / 



V^T^^TT^ P ^^E P1 + E p2 -E k -E/ 

1 d 3 l' 

W (^ow (5 - 3) 

We define the parton correlation function (PCF) f(xi ) as [12], 

f(xi, X2] x[, x' 2 )S(xi + x 2 — x[ — x' 2 )dxidx' 1 dx2dx' 2 
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E P1 + E V2 JE v > x + Ep> 2 



M p ^ PlP2 x [M p ^, p > iP > 2 x 



_Z_yp l^pj^ J^ p2 J^A 

1 1 1 1 „ d 3 A; x 



2E P1 2E P2 2E p , 2E P , (2tt) 3 2£ x ' 
while the parton recombination function Pj 2-2 -* is defined by 

dl 



(5-4) 



D (2-2)/ / / \i ™1 

n ; (xi,x 2 ,a; 1 ,a; 2 ,X3,X4)rfa;4^4- 

.Ei. 



-M PlP2 ^ k i> [Mp^p'^w 



E^ 2 JE^ 2 V ^ ™^ ^E P1+ E P2 -E k - Ei 



1 , 2 d 3 /' 



W (2vr) 3 2P/ (5 ' 5) 

We shall discuss the PCF and the parton recombination function in sections 6 and 7, 
respectively. Therefore, we have 

d<f)(x B , Q 2 ) 
d\nQ 2 

= (~^q) 2 J f(xi,x 2 ,x' 1 ,x' 2 ,Q 2 )8(x 1 + x 2 -x[ - x' 2 )Pf~ 2) {x u x 2 , x[,x' 2 ,x 3 , x 4 ) 

S(x 3 — xb)S(xi + x 2 — x 3 — x \)dx\dx 2 dx' x dx' 2 dx 3 dx 4, (5.6) 

We shall discuss the PCF and the parton recombination function in sections 6 and 7, 
respectively. Therefore, we have 

dq(x B ,Q 2 ) 

ding 2 

1 t 

= (~ftq) 2 J f(x 1 ,x 2 ,x' 1 ,x' 2 ,Q 2 )5(x 1 +x 2 -x[ - x 2 )Pi 2 ~ 2) (x l ,x 2 ,x' 1 ,x 2 ,x 3 ,x A ) 

S(x 3 — xb)S(xi + x 2 — x 3 — x ^)dx\dx 2 dx' x dx' 2 dx 3 dx ^ (5.5) 

where we have inserted a factor 1 — J 8{x\ + x 2 — x 3 — x±)dx 3 . This is the evolution 
equation from fig. 5a. 

Similarly, using the factorization in DIS [9,10], the contribution of fig. 5c (virtual 
diagram) is 

d(f>(x B ) = \-prd°{l*P -> k'X) 
Z O0 

1 \J Epi + E P2 \J E v\ + E P ' 2 
= 2 ^ Mp^ piX [M p ^p, iP2 p 2 x\ 



1 1 1 1 d 3 k x 



E p — E Pl — E x ^ E p — Ep^ — Ep 2 — E P2 — E x ^ 2E Pl 2E P2 2E p > i 2E p > 2 (27r) 3 2Px 
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H(-f*PiP2 -> l*PiP 2 )- (5-7) 

Now the PCF is defined by 

f(xi, x 2 , x[, x' 2 )dxidx[dx 2 dx 2 



M p ^ piX [M p -+ p > iP > 2P2 x\ 



E p 

1 , 1 



E p — E Pl — Ex E„ — E n [ — E n i_ — E nr , — E) 



1 1 1 \ ^ d A k 



2£ Pl 2£ P2 2£ p > Y (2n) 3 2E x 
However, we have the condition, 

f(xi, x 2 , x' 1 ,x' 2 ) = f(xi, x 2 ; x' 1: x' 2 ), (5.9) 

since the PCFs with cuts at different places are the same on the light-cone (fig. 6) 
[13,14]. 

The hard part is given by 

H(i*pip 2 -> 7V1P2) 

1 1 1 d 3 h' 

= (^) 2 ^^\M rP1 ^ 2 (^mp,+p 1 -k> 



1 1 d 3 l' 



^E pi +Ep 2 ^Ep,+Ep, 2 2E P1 (2tt 



! -n-zr-nr^-B eH*(^)*(^)* 



-Epi + -Ep 2 £V -Efc + Ej/ — Ep^ — E p ' 2 2Ek 2E\> 

' ' ' ' |M 7 * Pl ^(27r) 4 5 4 (p 7 + Pl -A;') ^ 



K R' C^Ep^ 10^^ > ^ ^ y (27r) 3 2^ 
£/ 1 

i= + ^ + ^ "^[M^n- - ^ - Efc - - r 

(5.10) 



v 2£ fc 7 (2vr) 3 2£/ 

where the factor of \ in (5.7) is needed for the cancellation of IR-divergences and 
momentum conservation as we shall discuss shortly. One can re-understand this factor as 
follows: only half of the probe-vertex connects with the partonic matrix in figs. 5c-f as 
well as in figs. 4b-c, and the square root of the parton density accepts the contributions 
of the partonic processes through a parton line. That is, 

y/qdy/q = 1 dq 

d\nQ 2 2d\nQ 2 ' K ' ; 
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Since, for the given initial partons we have, 

E M piP2 _ kll [M p[p ^ kl ,}* > 0, (5.12) 

s,t,u 

we can conclude that the negative sign in eq. (5.10) arises from 

' 1 ' f. (5.13) 



E Pl + E P2 — E k — Ei> E k + Ei' — E^ — E p > 2 E Pl + E P2 — E k — E, 
In consequence, we have 

d<j){x B ,Q 2 ) 1 1 



J f(x 1 ,x 2 ; x[, x' 2 , Q 2 )5(x 1 +x 2 - x[ - x' 2 ) 



dlnQ 2 2 y RQ' 

p(2-2)/ / / s 

5(a;i — xb)5(xi + x 2 — ^3 — x ^)dx\dx 2 dx\dx' 2 dx 3 dx ^. (5-14) 

Comparing eqs. (5.14) with (5.6), we see that the same recombination function appears 
in both cases. We can calculate the contributions of figs. 5a-f in TOPT using the same 
method, and finally obtain the hard contribution as 



d(f)(x B , Q 2 



Jq) 2 J f( x i^ x 2;x[,x' 2 ,Q 2 )Pi 2 2) {x 1 ,x 2 ,x' 1 ,x' 2 ,x 3 ,x A ) 



dhiQ 2 y RQ 
5(xi + x 2 - x[ - x 2 )[5(x 3 - x B ) + 5(x 4 - x B ) - ^5(xi - x B ) 

~\ S ( x i ~ Xb ) ~ \^ X2 ~ Xb ^> ~ \ S ( X 2 - x b)\ 

S(-Xi H — x[ H — x 2 H — x' 2 — x 3 — x A )dxidx 2 dx' l dx' 2 dx 3 dx 4. (5.15) 
2 2 2 2 

Obviously, eq. (5.14) contains the momentum conservation condition: 

d [ x B (p(x B ,Q 2 )dx B 

J ° Jln q. = °- < 5 ' 16 > 

This completes the discussion of the first type of diagrams. As the next step, we discuss 
the interference terms, iW^*_; fc /j/_ x -[M^,,_; fc /j/_ x -]*, shown in fig. 7. Proceeding similarly, we 
obtain the contributions from the interference processes in fig. 7 to be 



d(j)(x B , Q 2 



= 2P inter (— V q) 2 J f(x 1 ;x 2 ,x[,x' 2 ,Q 2 )Pi 1 3 \x 1 ,x 2 ,x[,x' 2 ,x 3 ,x A ) 



dlnQ 2 ul ^' y RQ 

5(x 1 - x 2 - x[ - x' 2 )[5(x 3 - x B ) + 5(x 4 - x B ) - ^S(x 1 - x B ) 
-\K x 'i - x b) - ^S(x 2 - x B ) - ^S(x' 2 - x B )} 
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S(x\ — X3 — x^dx\dx 2 dx' x dx' 2 dxzdx^ (5-17) 

where 



E pi + E p2 J E p i^ + E p ' 2 

f(x 1 ;x 2 ,x' 1 ,x' 2 ,Q 2 ) = - M p ^ plX [M p 

hi. 



*p\p' 2 P2^\ 



V 



1 1 1 1 „ d 3 k x 



-( = ) — — - — - — —IF 



E p — E Pl — E x > E p — Ep> x — E p > 2 — E P2 — E x 2E Pl 2E P2 2E p > i 2E p > 2 x (2%) 3 2E X 

H( 1 * PlP2 - 7 yy 2 ). (5.18) 



and 



J/2 

P 4 (1_3) (xi, x 2 , x[,x 2 , X 3 , X^dXA-^- 
Eh 



\J E Pl + E P2 \J E p[ + E p' 2 



-.M Pl ^ kV [M p ' iP ' 2P2 ^ k ii] 



1 1 1 d 3 I' 

<^) 2 jALr- (5-19) 



E k + E [ — E Pl E k + E v — E P [ — E P ' 2 — E p 2 (2tt) 3 2Eii 
In eq. (5.17), Pinter = or 1 implies that the interference processes are inhibited or 
exhibited, respectively. Now an interesting observation is that P 4 2-2 ^ and P 4 X ^ are really 
similar except for a simple coefficient and the different variable range. In fact, from fig. 
8 we find 

p (2-2) = _^p 4 (i-3) > ( 5 20 ) 

where £ is defined below and the negative sign occurs because (2xip) — l\j (2xi>p)) 
changes its sign from x\ < xy to x\ > xy in fig. 8. The contributions of the vertices A 
and B have same form, since the momenta of the partons a and b are [x a p, 0, x a p] and 
[xbp, 0, Xbp] , respectively. The factor £ arises from the following symmetry: if both the 
final partons are gluons or quarks (we do not distinguish quarks and antiquarks in this 
work), the corresponding virtual diagrams in figs. 7c-f are symmetric under the exchange 

(2—2) 

of these two partons. However, this symmetry will be lost if we use P 4 to replace 
P^ 1 ' 3 ^. In this case, £ = 1/2, otherwise, £ = 1. 

It seems that there are different energy deficits in going from P 4 X ^ to P 4 2 ^ in the 
ATOPT: l/(E k + E v - E pi ) in P 4 (1 ~ 3) and l/(E k + E v - E pi - E P2 ) in P 4 (2 ~ 2) ; l/(E k + 
E v - P p ; - E v i 2 - E V2 ) in P 4 (1_3) and l/{E k + E v - P p ; - E^) in P 4 (2_2) . However, they 
are really the same factor, arising from the term, (l\/(2x^p) + t]_/{2x4p))~ 1 . 

Therefore, one can replace P 4 X ^ by P 4 2 ^ and reconstruct eq. (5.17), where the terms 
f(xi] x[, x' 2 , x 2 , Q 2 ) and 5(x 1 — x[ — x' 2 — x 2 ) should be replaced by f( y x 1 ,x 2 ;x' 1 ,x' 2 ,Q 2 ) 
and S(xi + x 2 — x[ — x' 2 ), respectively, since Xi are the scaling variables. 

The final results from fig. 7 are 

d<l>{x B ,Q 2 ) 
d\nQ 2 
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1 r 

= 2(— ) 2 P inter J f{x l ,x 2 ;x' 1 ,x 2 ,Q 2 )Pi 2 ~ 2) {x l ,x 2 ,x' 1 ,x 2 ,x 3 ,x A ) 
8(x 1 + x 2 - x[ - x' 2 )[-8(x 3 - x B ) -8(x A - x B ) + ^S(x 1 - x B ) 

+^( x 'i ~ x b) + \K X 2 - x B ) + ^8(x' 2 - X B )} 

1111 

8(-Xi + -x[ + —X2 + -x' 2 — x 3 — xi)dx\dx2dx' x dx 2 dx 3 dx{. (5-21) 

Obviously, the momentum conservation condition is also satisfied in eq. (5.21), 

d ( x B (f)(x B ,Q 2 )dx B 

- ^v- = a < 5 - 22 ' 

6 Parton correlation functions 

In general, the parton density is a concept that is only defined at the twist- 2 level; it 
can be expressed in terms of the product of the initial and final hadronic wave functions 
with the same parton configuration. The parton correlation function f(xi, x 2 ; x[, x 2 ) is 
a generalization of the parton density beyond the leading twist. It has not yet been ex- 
perimentally observed. In this section, therefore, we shall try to construct the connection 
between the parton correlation function and the parton density. 

For example, consider the correlation function for the case when x\ = x[ and X2 = x' 2 
in eq.(5.8); this is given by 

f(xi,x 2 ;x 1 ,x 2 ) 

E Pl E P2 2 1 2 1 2 1 \ 2 1 r d 3 kx 

E p |M ™> 2X| [ Epi + E P2 +E X -E P } [ 2E P1 } [ 2E P2 } [ } (2nf2E x 

= p(x 1 ,x 2 ), (6.1) 

and is the number density of two partons, i.e., the probability of simultaneously finding 
two partons carrying x\ and X2 fractions of the proton momentum respectively. In the 
quantum mechanics approximation, we can use wave functions to represent p{x\,x 2 ) as 

p(x 1 ,x 2 ) = ip(x 1 ,x 2 )ip*{x 1 ,X2), (6.2) 

where ip(xi,X2) is the wave function of two partons in the proton. Similarly, we 
express the parton correlation function as the product of the initial and final hadron wave 
functions with different parton momentum, that is 

f(x u x 2 ; x[, x' 2 ) = if>(xi,x 2 )il>*(x' 1 ,aQ. (6.3) 

We define 

J(x u x 2 ;x' 1 ,x' 2 ) =ip(^ 1 ,x' 2 )ip*(x 1 ,x 2 ) = K 2 f{x x ,x 2 \x' x ,x 2 ). (6.4) 
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Therefore, we have 



f(x 1 ,x 2 ;x[,x' 2 ) = K~ 1 y/f(x 1 ,x 2 ; x[, x' 2 )J(xi, x 2 \ x[, x' 2 ) 

= K~ 1 yJtf(x 1 ,X2)^*(Xi,X2)tf(x[,X2)tf*(x[,X 2 ) 
= K~ 1 y/p(xi,X 2 )p(x' 1 ,X 2 ). 

In general, the two-parton number density can be approximated by 



(6.5) 



(6.6) 



p(x 1 ,x 2 ) = q a (xi)q b (x 2 )R 2 ab (x 1 ,x 2 ), 

where R 2 b (xi, x 2 ) is the momentum correlation of two initial partons; q a {x\) and qb(x 2 ) 
are the parton number densities. 

In order to estimate the value of k in eq. (6.5), we consider the process shown in fig. 
9. The time reversal invariance requires that 



da(j*p — > 7*p) = da(j*p <— 7*p), 



(6.7) 



or 



ip(x u x 2 )AB*^*(x[, x' 2 ) = ^K, x' 2 )BA*i/j*(x u x 2 ), (6.8) 
where A and B are the contributions of the hard parts in fig. 9. Therefore, 



k = 



/(ari,x 2 ;xi,X2) 
\ f(x 1 ,x 2 ;x' l ,x' 2 ) 



i AB* 
~BA* 



(6.9) 



Since k is expressible in terms of hard parts, eq. (6.9) indicates that k is calculable 
within PQCD. 



7 New evolution equations 

We now apply the method, used to describe scalar partons in Section 5, to the realistic 
case of partons (quarks and gluons) interacting within QCD. In consequence, we have 
following new evolution equations with twist-4 for GG — > qq and GG — > GG respectively: 
(a) GG — > qq. The contribution to the evolution equation for gluons is, 

dG(x B ,Q 2 ) 
dlnQ 2 

= (^) 2 / y/Gfa, Q 2 )G(x 2 , Q*)G( Xl + A, Q 2 )G(x 2 - A, Q*) 
Rgg(xi,x 2 )R gg ( Xi + A, x 2 - A) 

J2^T 1P GG^qq( X l^ X 2,X 3 ,X 4 , A)S(x 1 + X 2 - X 3 - X 4 ) 
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11 11 

[—-S(xi — Xb) — 77^(^1 + A — Xb) — 77^(^2 — Xb) — 77^(^2 — A — X B )]dXidx 2 dx 3 dXidA 

jL Zj Zi Zi 

+2{ ^Q )2Pmter J ^ G{XU Q ^ G ^ Q 2 )°^ + A > Q 2 )G{x2 - A, Q 2 ) 
Rgg(xi, x 2 )Rgg{x 1 + A,x 2 - A) 
^k^P^^^i, £2,^3,^4, A)5(xi +x 2 -x 3 -x 4 ) 

i 

11 1 1 1 

-[-5(x 1 -x B ) + -5(x 1 + A-x B ) + -5(x 2 -x B ) + -5(x 2 -A-x B )]dx 1 dx 2 dx 3 dx 4 dA, (7.1) 

where the factor | in the last factor arises from symmetry considerations, just as in 
eq. (5.20). However, this symmetry will be broken due to the cut in real diagrams in the 
corresponding equation for quarks: 

dq(x B ,Q 2 ) 
d\nQ 2 

= (^) 2 / y/G(*i, Q 2 )G(x 2 , Q 2 )G( Xl + A, Q 2 )G(x 2 - A, Q 2 ) 

R GG (x 1 ,x 2 )R GG (xi + A, x 2 - A) 
k~ x P GG ^ q y(xi, x 2 , x 3 , X4, A)5(xi+x 2 — x 3 — X4)[S(x 3 — x B )+5(x 4: — XB)]dxidx 2 dx 3 dx±dA 



+2(j^) 2 P tnter f y/G(x u Q 2 )G(x 2 , Q 2 )G( Xl + A, Q 2 )G(x 2 - A, Q 2 ) 

R GG (x 1 ,x 2 )R GG (xi + A, x 2 - A) 
^ j K~ 1 P l GG _ ¥m (x 1 ,x 2 ,x z ,x 4 ,A)8(x 1 +x 2 -x 3 -x 4 ) 

i 

[— 5(x 3 — x B ) — 5(x4 — x B )]dxidx 2 dx 3 dx4dA. (7.2) 

(b) GG — > GG. The contribution to the evolution equation for gluons is, 

dG(x B ,Q 2 ) 
dlnQ 2 

J y/G( Xl , Q 2 )G(x 2 , Q 2 )G( Xl + A, Q 2 )G(x 2 - A, Q 2 ) 

R GG (x 1 ,x 2 )R GG (x 1 + A, x 2 - A) 
^2Ki 1 PQ G ^ GG (xi,x 2 ,x 3 ,X4,A)S(xi + x 2 -x 3 -x 4 ) 

i 

[5(x 3 -x B ) + 5(x4-x B )-^5(x 1 -x B )-^5(xi + A-x B )-^5(x 2 -x B )-^5(x 2 -A-x B )] 

dx\dx 2 dx 3 dx4 dA 

+2(^) 2 P tnter J y/G(x u Q 2 )G(x 2 , Q 2 )G( Xl + A, Q 2 )G(x 2 - A, Q 2 ) 
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R gg (x 1 ,x 2 )Rgg(x 1 + A, x 2 - A) 

^2 ^ P GG^Gg( X 1^ X 2, x 3, x 4, A)5(x 1 + X 2 ~X 3 - X A ) 

i 

11 11 

[-5(x 3 -xb)-S(x4-xb) + -S(x 1 -xb) + -S(x 1 + A+xb) + ^S(x2-xb) + ^S(x 2 -A-xb)] 

dx\dx 2 dx 3 dx 4 dA. (7.3) 
We discuss the case when Pinter = 1 and separately. 

A: Pinter = 1- We can take Pinter = 1 if there is no reason to forbid three-parton 
recombination in the interference terms in nucleon. Thus, we have 

dG(x B ,Q 2 ) 
dlnQ 2 

= (t^) 2 / y/G( Xl , Q 2 )G(x 2 , Q 2 )G( Xl + A, Q 2 )G(x 2 - A, (?) 

R G g(xi,x 2 )R G g(xi + A, x 2 - A) 

J2 K i 1P GG^GG( X l^ X 2, x 3,X4,A)S(x 1 + X 2 - X 3 - X 4 ) 

i 

[S(x 3 — xb) + 5{xi — x B )}dxidx 2 dx 3 dx 4 dA 
-2(^r) 2 / y/G( Xl , Q 2 )G(x 2 , Q 2 )G( Xl + A, Q*)G(x 2 - A, Q 2 ) 

tiX^) Jx\>x B 

Rgg(xi,x 2 )R gg (xi + A, x 2 - A) 

X] K 7 lp GG^GG( X l, x 2,%3, X 4 , A)5(x 1 + X 2 - X 3 - X A ) 



and 



[S(x 3 — x B ) + 5{x 4 — XB)]dxidx 2 dx 3 dxidA. (7.4) 

dq(x B ,Q 2 ) 
d\nQ 2 

■ (t^) 2 I y/G(x lt Q 2 )G(x 2 , Q 2 )G( Xl + A, Q 2 )G{x 2 - A, Q 2 ) 

rLLJ J (xi+x2)>x B 

Rgg(xi,x 2 )R G g(xi + A, x 2 - A) 
J2k~ 1 P gg ^(x 1 ,x 2 ,x 3 ,x 4 , A)S(x 1 +x 2 -x 3 -x 4 ) 

i 

[S(x 3 — x B ) + 5(x 4 — x B )]dx\dx 2 dx 3 dx 4 dA 

-2(^r) 2 / ^(xx, g 2 )G(x 2 , Q 2 )G( Xl + A, g 2 )G(x 2 - A, Q 2 ) 

Rgg(xi,x 2 )R gg (xi + A, x 2 - A) 
^K~ 1 P l GG ^ m (x 1 ,x 2 ,x 3 ,x 4 ,A)8(x 1 +x 2 -x 3 -x 4 ) 

i 
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[8(x 3 — xb) + 5{x 4 — XBf\dx\dx 2 dx 3 dx 4 dA. (7-5) 

Equations (7.4) and (7.5) predict that the shadowing effect in quark distributions is 
stronger than that in the gluon distribution, since there are two shadowing sources for 
quarks but only one shadowing source for gluons. 

B: Pmter — 0. This means that the interference terms are forbidden. An example of 
such a case is the radiation recombination in a nucleus. We consider the recombination 
of partons which originate from different nucleons in a nucleus. A single parton can not 
escape from the confinement region of a nucleon, unless it forms a colour-single cluster 
with other partons. We define the probability of a parton leaks out from the confined 
volume as w. Thus, 



and 



M {2 - 2 \x u x 2 ^x 3 ,x 4 )[M^ 2 ' 2 \x 1 ,x 2 ^x 3 ,x 4 )}* oc w, (7.6) 

M {2 ~ 1) (x 1 -> x 3 ,x 4 )[M (2 - 3 \x 2 ,x[,x' 2 -> x 3 ,x 4 )\* 

+M {2 - 3) (x 2 ,x[,x 2 -> x^x A )[M {2 - l \x l ^2:3,0:4)]* ocw 3 , (7.7) 

We can neglect the interference processes (7.7), because of the confinement condition 
w < 1. In this case, we have another face of the evolution equation: 

dG(x B ,Q 2 ) 
d\nQ 2 

= (^) 2 / y/G(x lt Q 2 )G(x 2 , Q*)G( Xl + A, Q*)G(x 2 - A, Q*) 

R GG (x 1 ,x 2 )R GG (x 1 + A, x 2 - A) 
^k~ 1 P gg _ gg (x 1 ,x 2 ,x 3 ,x 4 , A)5(x 1 +x 2 -x 3 -x 4 ) 

i 

[S(x 3 -xb) + S(x 4 -xb)-^S(x 1 -xb)--S(x 1 + A-xb)-^S(x2-xb)--S(x 2 -A-xb)} 

dx\dx 2 dx 3 dx 4 dA 
~(^q) 2 J y/G(x u Q 2 )G(x 2 , Q 2 )G{ Xl + A, Qi)G(x 2 - A, Q 2 ) 
Rgg(x 1 ,x 2 )R gg (x 1 + A, x 2 - A) 

Yl K i 1P GG^qq( X l^ X 2,X3,X 4: A)S(x 1 + X 2 - X 3 - X 4 ) 

i 

[^S(xi - x B ) + \ 5 ( x i + A-x b ) + ^5{x 2 - x B ) + ^{x 2 - A - x B )} 

dx\dx 2 dx 3 dx 4 dA. (7.8) 

and 
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dq(x B ,Q 2 ) 
d\nQ 2 

= (^) 2 / y/G(xu Q 2 )G(x 2 , Q 2 )G(x l + A, Q 2 )G(x 2 - A, Q 2 ) 

Rgg(xi,x 2 )R gg (x 1 + A, x 2 - A) 
E K i lp GG^qq( x i , x 2,x 3 ,x 4 , A)5(x 1 +x 2 -x 3 -x 4 )[5 (x 3 -x B ) +5(x 4 -x B )]dx 1 dx 2 dx 3 dx 4 dA 



(7.9) 

Now the sign of the right-hand side of eq. (7.9) is positive. This means that the 
shadowing in quark distribution is weaker than that in gluon distribution. 

In principle, we can calculate the parton recombination functions at order 0(a 2 ) for 
every parton flavors in the whole x region. However, in the majority of cases, the parton 
recombination is happened in the gluons with small x. For simplicity, we only consider 
the case where all partons are gluons with small x value in this work. We will discuss the 
recombination of partons in a general x range elsewhere. In this approximation, we can 
use the results of Mueller and Qiu in the calculations of the real process of fig. la in ref. 
[2]. Thus, the contributions of figs. 6a,b to P 4 are from t- and u-channels and as well 
as their interference-terms. One obtains 

i 1 

where we assume that x\ = x[ = x 2 = x' 2 . Eq. (7.10) evaluates to 

2 J E p( GG ~ GG - } 5(x 3 — xb)5(2xi — x 3 — x 4 )dxidx 3 dx 4 

i 

4tt 3 a s C A 2 r dxi 

" N 2 -V 7T ' J Xl ' 1 ' 

which is the same as the result of ref. [2]. If we take the gluon correlation function [2] 
to be 

^G«(x 1> xi,x 2 y 2> Q 2 ) = ^G(x 1 ,Q 2 )G(x 1 ,Q 2 ), (7.12) 

we obtain the following simplified evolution equations arising from GG — > GG in the 
small x region: 

dx B G(x B , Q 2 ) _ 81,a s )2 [x B G(x B , Q 2 )] 2 



d\nQ 2 \Q y RQ J x B 

16 RQ Jx B /2 x{ 



, p 81 a s 2 [x B G(x B ,Q 2 )] 2 

+-T inter n „ V D< n / 

16 RQ x B 
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Note that the modifications of fig. 5a-b are related to G(x B ,Q 2 ) in eq. (7.13) but 
not to x b G(xb, Q 2 ) according to eq. (5.1) in our work. However, these real diagrams fig. 
5a-b (or fig. la) in ref. [2] are regarded as the modifications in xbG(xb,Q 2 )- Therefore, 
the equation (7.13) is different with the GLR equation in the dependencies of x\ and x B . 

In consequence, the new evolution equations if Pinter = 1 are 

dx B G(x B ,Q 2 ) _ 81a, f 1 / 2 [x 1 G(x 1 ,Q 2 )} 2 
dlnQ 2 16 W XB Lb/2 x\ 1 

81 ' fts %/ 1/2 ^, (7.14) 



8 RQ Jx B x 
Here, we have an extra conservation 



Jo 1 dx B x B G(x B , Q 2 ) _ 

2 = 0> (7.15) 



rl rl/2 rl/2 rl/2 

/ dxs / f(x\)dxi — 2 / dx B / f{xi)dx\ = 0. (7-16) 

JO Jx R /2 JO Jx H 



dlnQ 

since for any function f(xi) we have 

-1 r l/2 rl/2 rl/2 

f(xi)dxi - 2 / dx B / 

C B /2 JO Jx B 

On the other hand, the new equation has following different form if Pinter = 

dx B G(x B ,Q 2 ) 
dlnQ 2 

81, a, ^ 2 [x B G(x B ,Q 2 )] 2 81^ f 1 ' 2 [xMx^Q^l 
~ 16 W ' x B " + 1Q [ RQ } XB Jx B/ 2 ' x\ ' dXl ' U '> 

8 Discussions and conclusions 

The following interesting components of the new evolution equation derived in this 
paper are highlighted: 

1. Through the derivations of sections 4 and 5, it seems there is an interesting "cutting 
rule" in DIS: The contributions of the cut diagrams in the sum (3.4) have the 
identical integral kernel with only the following different factors R: 

R=(±)x(l,^)x5(xp-x B ). (8.1) 

The various terms appearing in the cutting rule (8.1) can be described in terms of 
the general structure of the cut diagrams G 7 (iV) in TOPT: 



g ^ n )= n ~ — rp f t? n 



left-vertices ^a'ei+1 ^a' J2aei right-vertices ^bej J2b'ej+1 ^b' 
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n o) n^n| n 7=^, <") 

vertices— k / states— c c loops— d V / final— states— e V / e 

where iVc is an overall numerator-and-symmetry factor and is independent of the 
cut 7 [6]; i and i + 1 (or j and j + 1) are the time-ordered lines on the left- (or 
right-) vertices; 6(52/ Xf — 0) is the conservation of longitudinal momentum at the 
vertex. 



(a) The sign in the first factor of (8.1) is determined by the energy deficits in (8.2). 
For example, if a vertex pass through the cut line, the corresponding energy 
deficit will change its sign since 



1 



(8.3) 



52a'€i+l E a ' 52aei 52a£i E a 52a'€i+l ^a' 

as we have for example in (4.11). 

(b) The second factor takes a value of 1/2 if the probe- vertex inserts in the initial 
line as shown in (5.11). 

(c) 8{xp — xb) is the direct result of the sum (3.4). 

(d) When the cut line moves its position, the contributions of the final states in 
(8.2) will change the momentum-symbols, but don't change the structure of 
the intermediate state according to the sum-condition I in (3.4). We also note 
that a virtual parton line has 4-dimensional integral and a real parton only a 3- 
dimensional integral in the covariant perturbation theory; however, in TOPT, 
since the virtual partons are replaced by the effective real partons, the above 
mentioned differences are contained in the energy deficits. In particular, when 
the cut line moves to cutting line a c from cutting a loop d-d! in the process of 
c — > d'd — > e, we have similar expressions in (8.2) due to 



d°k d , cfk rf d 6 K d 6 k d 1 1 



2E C {2iif2E dl {2nf2E d (2tt) 3 2£ c (2tt) 3 2E d 2E d , ' 



(8.4) 



Thus, as the cut line moves from cutting a loop to uncutting a loop, the integral 
kernel has the same form. The difference lies only in the cutting positions. 

(e) The cut line can also cut the nonperturbative matrix elements with multi- 
initial partons. The reason is that the initial parton line on the light-cone can 
be moved from one-side of the cut to another side (see fig. 6) [13, 14]. Thus, 
we can use the same correlation function to represent the different hadronic 
parts and keep the same integral kernel. Because of this important property, 
the sum (3.4) shall include more complex cut diagrams when we study parton 
fusion or recombination. 

(f) Finally, the matrix can been factorized to obtain the probability explanation; 
the reason is that we used the bare probe-vertex approximation and the co- 
efficient C q in (2.2) (or in (5.2)) is canceled in the calculations at this 
approximation. 
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2. In the DGLAP equation, there are IR-divergences when a final-state gluon becomes 
soft. These IR-divergences can be taken care of by the sum of figs. 4 according to 
eqs. (3.4) and (8.1). The divergences in (7.11) are canceled due to the symmetry 
at small x approximation [2]. However, a soft initial parton also may give rise to 
IR-divergences in the parton recombination process in a general case. We now show 
that such IR-divergences can be canceled by using the same method in the DGLAP 
equation. 

For example, take x 2 = in fig. 10; this implies p 2 = 0. Since the unpolarized 
structure functions only involve contributions from terms with even- twist, we have 
x' 2 = at x 2 = 0. The momentum of II and Ir are determined by the down- 
vertices since the energy is not conserved in the up- vertices in the TOPT, i.e., 

i 2 i 2 
h = [(x 2 -x 4 )P-^p,h, (x 2 -x A )P\ and l R = [(x 4 -x' 2 )P + ^jp, -l±, (x 4 -x' 2 )P\. 

Thus, x 4 = and x 1 = x[ = x 3 . According to (5.14), we can find that 

J f(x 1 = x B , x 2 = 0, x[ = x B , x 2 = 0, Q 2 ) 

[p 4 ( 2 - 2 )( Xl — X3? X2 — o ; x [ = x 3 , x 2 = 0, x 3 , x 4 = 0)5(x 3 - x B )dx 3 

— P^ 2 ~ 2 \xi, x 2 = 0, x[ = Xi,x' 2 — 0, x 3 — xi, X4 — 0)5(xi — x B )dxi 
2 

— P4 2 ~ 2 \xi = x[,x 2 = 0,x[,x 2 = 0,x 3 = x[,X4 = 0)5(x' 1 — x B )dx[] = 0. (8.5) 

Therefore, IR-divergences can be canceled point-by-point at the IR-pole. 

3. Obviously, the new evolution equations (7.4), (7.5), (7.8) and (7.9) are different 
from the GLR equation [1,2]. It is interest that the properties and structure of 
the simplified \ow-x form (7.14) is similar to a modified GLR equation that has 
been obtained earlier in ref. [15]. However, two equations really have a different 
theoretical basis. The GLR equation and its modified form are based upon the 
AGK cutting rules. They sum three kinds of diagrams: cutting two-ladders, one- 
ladder and zero-ladder, respectively. The first (see fig. la) is identical with our 
figs. 5a,b, however, the cut lines in the latter two figures (figs. lb,c) break parton 
recombination and these processes should be inhibited. 

In conclusion, parton recombination via QCD evolution equation is investigated using 
perturbative theory without the AGK cutting rules. The contributions from different cut 
and interference diagrams are summed and infrared safeness and momentum conservation 
are established. Time ordered perturbation theory is developed to establish the connec- 
tions among different cut diagrams. As a consequence, a new nonlinear evolution equation 
is derived on a different basis from the GLR equation. Furthermore, this new evolution 
equation is more detailed in structure than the GLR equation. 
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Figure Captions 

Fig. 1 Two expressions of the GLR equation based on the AGK cutting rules: (left) cut 
vertex theory and (right) time-odered perturbation theory; where the cuts through (a) 
two-, (b) one- and (c) zero-ladders, respectively. 

Fig. 2 The amplitudes (a) M^^l k , vx , (b) M^fl kll , x and (c) M^7^ k , l/X - The dark circles 

denote the PQCD interaction with the correlation of the initial partons. 

Fig. 3 Naive parton model of DIS. 

Fig. 4 The leading order splitting processes in DIS. 

Fig. 5 The diagrams contributing to the leading recombination order from \M^ t fl k , l , x \ 2 . 
Fig. 6 Identical hadronic parts in different cut graphs from refs. [13,14]. 
Fig. 7 The diagrams contributing to the leading recombination order from 

2 /i / f( 2 - 1 ) r/if (2 ~ 3) l* 

z ' lvl p~ 1 *->k'l'X[ lvl p~ 1 *^k'l'X\ ■ 

Fig. 8 A diagrammatic illustration of the relation between P] 2-2 " 1 and P^~^ ■ 
Fig. 9 A diagrammatic illustration of the time reversal invariance in 7^ — > 7p. 
Fig. 10 IR safeness in GG -> GG. 
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